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Abstract  ' 

Ribner’s  (1962)  dilatational  acoustic  theory  is  revisited.  A  rigorous  connection  between  the  fluctu¬ 
ating  dilatation  rate  and  the  acoustic  source  field  is  established;  this  vindicates  Ribner’s  heuristic 
contention  while  indicating  additional  acoustic  source  terms  in  his  dilatational  acoustic  theory.  It 
is  also  shown  that  Ribner’s  acoustic  source  term  is  quadrupole.  Interesting  consequences  of  the 
dilatational  point  of  view  are  indicated.  It  is  found  that  in  the  region  of  vortical  fluid  motion  the 
dilatation  scales  as  the  square  of  the  turbulent  Mach  number,  M^,  and  has  little  to  do  with  the 
acoustic  field;  its  time  rate  of  change,  however,  is  a  portion  of  the  sound  generation  mechanism. 
Away  from  the  vortical  region  the  fluid  dilatation  is  an  acoustic  variable  and  scales  as  Mf.  The 
mathematical  link  established  is  useful  to  interpreting  direct  numerical  simulation  of  aeroacoustical 
flows  in  which  the  dilatation  is  computed. 
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1.  Introduction 

In  the  nascent  field  of  direct  computation  of  aeroacoustical  flows,  the  dilatational  fields  are  some¬ 
times  investigated,  see  for  example,  Colonius,  Lele,  and  Moin  [1]  or  Mitchell,  Lele,  and  Moin  [2]. 
One  has,  in  general,  an  intuitive  appreciation  of  the  fact  that  sound  and  dilatation  are  linked.  This 
link  has  not  been  rigorously  confirmed  for  the  aerodynamic  source  problem.  Links  have  been  made 
to  fluctuating  Reynolds  stresses,  Lighthill  [3],  Phillips  [4],  the  fluctuating  vorticity,  Powell  [5],  Howe 
[6],  the  pseudo-pressure,  Ribner  [7],  the  fluctuating  triple  velocity  gradient,  Lilley  (1974)  (more  eas¬ 
ily  found  in  Goldstein  [10])  and  the  fluctuating  enthalpy  Howe  [6],  This  article  also  provides  and 
clarifies  some  mathematical  details  missing  in  [7]’s  heuristic  analysis.  In  addition  to  formally  es¬ 
tablishing  the  link  between  the  dilatation  rate  and  sound  source  some  of  Ribner’s  contentions  are 
vindicated.  Useful  DNS  consequences  of  the  dilatational  point  of  view  are  indicated.  As  the  dilata¬ 
tional  approach  to  an  aeroacoustic  theory  appears  to  have  generated  some  controversy  some  three 
decades  ago,  it  seems  necessary  to  make  clear  at  the  outset  that  the  object  of  this  article  is  not 
related  to  the  use  of  the  dilatational  theory  as  an  acoustic  analogy. 

In  the  present  era,  some  useful  items  related  to  Ribner’s  [7]  dilatational  theory  suggest  a  reap¬ 
praisal  of  Ribner’s  ideas.  Ribner  [7]  related  the  acoustic  source  to  the  second  time  derivative  of 
the  incompressible  pressure.  Ribner  [7]  gives  physical  arguments  relating  his  sound  source  to  the 
fluctuating  dilatation.  His  arguments  have  not,  however,  been  given  formal  mathematical  expres¬ 
sion  or  justification.  This  article  formally  establishes  Ribner’s  [7]  heuristic  assertion:  that  sound 
generation  is  related  to  the  time  rate  of  change  of  dilatation  The  derivation  of  this  mathematical 
connection  and  related  consequences  as  well  as  some  useful  aspects  of  the  dilatational  point  of  view 
are  the  subjects  of  this  article.  In  addition  to  exploring  and  clarifying  Ribner’s  [7]  ideas  some  useful 
and  relevant  consequences  for  DNS  and  LES  are  obtained.  In  focussing  on  Ribner’s  analysis,  the 
discussion  is  limited  to  compact  flows  with  negligible  nonisentropic  or  diffusive  effects.  Implicit  in 
the  development  is  the  limitation  to  the  class  of  flows  for  which  the  Lighthill  form  of  the  acoustic 
analogy  is  useful.  Many  of  the  effects  accounted  for  in  more  advanced  acoustic  theories,  such  as 
those  of  Lilley  or  Phillips,  are  not  addressed  in  this  article. 

In  summary:  1)  a  formal  link  between  the  sound  source  and  the  dilatation  is  derived;  2)  Ribner’s 
source  term  is  shown  to  be  quadrupole;  3)  additional  source  terms  in  Ribner’s  dilatational  theory 
are  determined;  and,  4)  the  potential  of  the  time  rate  of  change  of  the  dilatation  as  a  diagnostic 
tool  is  indicated. 

^For  this  article  the  quantity  Ui,i  will  he  called  the  dilatation  (it  does  in  fact  have  units  of  a  rate). 
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2.  Some  preliminaries  and  the  Ribner  analysis 

A  few  introductory  comments  and  an  informal  outline  of  some  useful  physical  ideas,  following 
Ribner  [7],  are  summarized.  In  this  section  Ribner’s  [7]  mathematical  development  is  followed;  this 
is  in  order  to  highlight  its  differences  with  a  more  mathematical  analysis. 

For  the  present  class  of  flows  the  term  (p  —  c^p),tt,  following  [7],  is  neglected  and  the  pressure 
fluctuations  in  a  fluid  satisfy  the  following  Lighthill  equation, 

-  P,3j  =  {p-UiUj),ij  (1) 

where  p  represents  the  deviations  of  the  fluid  pressure  from  a  static  reference  value.  In  medium  of 
uniform  mean  density  pso,  it  is  customary  to  approximate  pUiUj  =  since  higher  order  terms 

scale  as  the  square  of  the  fluctuating  Mach  number.  Here  Vi  represents  the  solenoidal  velocity. 
Ribner  [7]  assumes  at  the  outset  the  validity  of  this  approximation.  Ribner  distinguishes  two 
pressures,  an  “acoustic”  pressure  which  propagates  and  a  “pseudo-pressure”  ®  associated  with 
the  convective  motions  of  the  fluid.  The  pressure  field  in  the  Lighthill  equation  includes  both 
acoustic  and  convective  pressures.  The  fluid  pressure,  following  [7j,  is  partitioned  into  convective 
and  propagating  parts  p-=  p^  +  Pa  where  ps  satisfies 

—  V  Ps  =  Pco(^Ab)?b  •  (^) 

The  subscript  s  is  used  to  denote  the  incompressible  or  pseudo-pressure  associated  with  the 
solenoidal  velocity  field.  The  acoustic  pressure,  pa,  following  Ribner  satisfies 

Pa-U  -  V  =  -  (3) 

the  sound  source  has  been  related  to  the  incompressible  pressure  which  satisfies  the  Poisson 
equation  whose  solution  is  given  by  convolution.  Ribner  [7]  then  goes  on  to  argue,  without  giving 
proof,  that  Ps  (the  overdot  signifies  a  time  derivative)  is  related  to  the  dilatation. 

The  following  section  establishes  a  rigorous  connection  of  Ribner’s  contention  of  a  relation  between 
the  pseudo-pressure  sound  source  and  the  dilatation.  Some  subtleties  regarding  Ribner’s  very  useful 
decomposition  p  =  p^  -f-  p^  are  described.  It  is  seen  that  the  decomposition  retains  terms  that  are 
of  the  same  order  in  as  terms  discarded  in  the  truncation  puiUj  =  pa^ViVj. 

3.  A  perturbation  analysis  leading  to  Ribner’s  acoustic  analogy 

The  following  analysis  is  understood  to  address  high  Reynolds  number,  weakly  compressible  (low 
fluctuating  Mach  number),  turbulent  compact  flows  with  constant  mean  density  the  class  of  flows 
®The  term  “pseudo-pressure"  coined  by  Blokhintsev  (1956)  as  quoted  in  [7]  is  adopted. 
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relevant  to  subsonic  aeroacoustics  for  which  the  energy  equation  implies  p  ~  p^.  A  low  fluctuating 
Mach  number  perturbation  analysis  of  the  compressible  Navier-Stokes  equations. 


Pit  +  M/tP./t  =  -  pVkik 

(4) 

pviit  +pVkVi,k  +Pii  =  0 

(5) 

p/poo  =  {P/Poop 

(6) 

Pitt  ~  Pijj  =  (p''^j®j)-y 

(7) 

is  required.  Perturbing  about  a  reference  state,  (Pco-Poo)-  the  nondimensional  forms  of  the  pressure 
and  density  are  taken  as  p  =  Pco{l  +  p'),  p  =  Pco(l  +  p')-  The  independent  variables  are  rescaled 
with  the  energy  containing  length  and  time  scales  of  the  turbulent  flow  field:  and  1.  Dropping 

primes  the  compressible  Navier-Stokes  equations  become 


P,i+V-kP,k  =  -(1 (8) 

(l  +  p)«i,(  -t- (l-|-p)UpUi,p  +e”^p,i  =  0  (9) 

P-7P  =  1/2  7(7-1)p^  (10) 

.  Pm  -  (■~'^P,j3  =  [(1  +  •  (11) 

A  wave  equation  related  to  the  last  equation,  in  nondimensional  units,  is 

=  [(1  +  p)tti«:j],.y  +  [7“V  -  P]m-  (12) 

In  these  equations  =  7«?/c^  and  =  ypc^/poo-  The  characteristic  fluctuating  velocity 

might  be  identified  with  the  energy  of  the  turbulent  fluctuations:  ~  |A:  =  1  <  UjUj  >.  The 


conventional  definition  of  the  Mach  number  (in  contradistinction  to  the  larger  fluctuating  Mach 
number  used  in  DNS  of  compressible  turbulence)  is  used.  Expansions  of  the  form  p  =  [  Pi  +e^P2+ 

...],  p  =  [  Pi  -t-  e^p2  -h  ...],  Vi  =  t’i  +  [  Wi  -1-  +  •••]-  are  chosen.  Substituting  the  expansions 

into  the  equations  produces,  to  the  lowest  order,  the  incompressible  form  of  the  equations. 


+  Plri 

=  0 

(13) 

t-’i-i 

=  0 

(14) 

V"pi 

=  - 

(15) 

Pi 

=  7Pi- 

(16) 

Note  that  variables  are  now  dimensionless.  The  lowest  order  form  of  the  equations  of  state  indicates 
that  the  pressure  is  the  solenoidal  or  pseudo-pressure  ■*.  Additional  details  can  be  found  in  Ristorcelli 

[13L _ 

^In  dimensional  form  this  is  statement  that  pi  =  c^paj pao- 
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The  inner  form  of  the  wave  equation  is,  to  order  Mf, 

-V^P2  =  [(l+Pi)«iWi  -  ViVjUj  -y~^Putt-  (17) 

Note  that  for  brevity  and  clarity,  the  quantity  UiUj  has  not  been  expressed  in  the  pertubation 
variables.  The  term  in  square  brackets  is  more  physically  meaningful  if  left  as  is:  it  is  the  difference 
between  the  compressible  and  incompressible  fluctuating  Reynolds  stresses.  The  continuity  equa¬ 
tion,  at  next  order,  is  pi,t  +  —y>h<k-  Using  the  equation  of  state  pi  =  'ypi,  the  first  order 

continuity  equation  yields  a  diagnostic  relationship  for  the  leading  order  dilatation,  d  =  W/;,/;, 

-yd  =  pi,(  +VkPi,k.  (18) 

Ribner’s  [7]  interpretation  of  the  second  time  derivative  of  the  pseudo-pressure  as  dilatation  rate  is 
now  formally  established:  pi  =  —yd  —  {vkPk)<t  ^-nd  P2  satisfies,  to  0{M^), 

-V'^p2-  d+  [(l  +  pi)ttiUj  -  WiVj],,:,-  -|-7”^(nfcpi),«.  (19) 

Lighthill’s  equation,  C^p.u  —  P^jj  —  (p'Ui'Uj),ij,  can  be  reconstituted  in  dimensional  terms  (with¬ 
out  implied  length  and  time  scales).  Using  Ribner’s  decomposition  p  =  Ps  +  Pa  where  the  dwnen- 
siona/ quantities  [ps,p.j]  corresponding  to  the  nondimensional  [pi.ps]  produces,  to 

c7jpa,U  -  VV.»=  7C«Poorf+  Poo[(l+p')«i«i  -  (20) 

The  divergence  of  the  pressure  flux  is  an  octupole  (see  Appendix)  and  of  higher  order  in  than 
the  dilatation  rate.  The  wave  equation  in  Ribner’s  dilatational  acoustic  theory  is,  consequently, 

cTjpa.tt  -  =  7ciV.<,rf  +  pc»[(l+p')«.“j  -  (21) 

The  dilatation  is  seen  to  be  one  portion  of  the  acoustic  source;  additional  quadrupoles  (in  the  square 
brackets)  arise  from  the  compressible  nature  of  the  Reynolds  stresses.  These  terms  do  not  appear 
in  the  Ribner  analysis  as  pUiUj  ~  pooViVj  was  assumed  at  outset  -  a  truncation  of  terms  of  the  same 
order  as  d.  (The  source  term  can  be  expanded  [(1 +p')«,-«j  —  UiUj]  =  ViWj-\-VjV)i+ p'viVj+h.o.t.) 

It  should  be  noted  out  that  Crow’s  [9]  analysis  of  Ribner’s  theory,  in  which  it  is  shown  that  Ribner’s 
source  term  /  Psdy  =  0,  would  have  been  different  if  the  retarded  time  had  been  included.  Retention 
of  the  full  retarded  time  would  have  given  a  non-zero  integral  of  Ribner’s  source  term  invalidating 
one  of  Crow’s  conclusion  regarding  Ribner’s  theory.  Ffowcs  Williams  [12]  indicates  that  ps  ~ 
and  the  integral  is  in  fact  weakly  (logarithmically)  divergent  requiring  very  careful  source  term 
modeling  and  integration,  if  it  is  to  be  used  as  an  acoustic  analogy.  In  contradistinction,  the  source 
term  in  the  Lighhill  analogy  scales  as  the  double  divergence  of  ViVj  ~ 
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While  the  expression  c^Pa.u  —  V^Pa  =  yc^Pxd  +  — •  identifies  a  new  quantity,  d,  it  does  not 
contain  any  new  physics.  The  quadrupole  nature  of  the  sound  source  field  is  retained  as  can  be  seen 
from  the  following  application  of  the  convolution  algebra;  The  solution  for  the  solenoidal-pressure 
is 


Ps{x)  =  [uiUj (*')], i-y  ®G'o(a;-  »')  =  (viVj)  =  [viVj)  ®  G o,ij  =  SiGolvij  •  (22) 

Go  is  the  three-dimensional  Greens  function  for  the  Laplacian,  Go  =  ~  *'1“^.  The  dilatation 

rate  is  related  to  the  second  time  derivative  of  the  incompressible  pressure  and  the  source  term  is 
re-expressed  as 

d  ~  Pj  ~  ®  Clo],ij  ■  (23) 


In  less  abstract  terms 


dx' 


|x- 


lij  — 


(24) 


The  source  term  is  a  double  divergence  and  the  quadrupole  nature  of  the  source  term  is  retained. 
It  appears  that  d  (and,  of  course,  Ribner’s  original  p,)  are  to  be  understood  as  a  quadrupole. 


In  summary:  In  this  revisit  to  Ribner’s  [7]  theory  1)  the  source  term  p^,  has  been  formally  linked 
to  the  dilatation;  2)  the  source  term  in  his  theory  [7]  has  been  shown  to  be  a  portion  of  the  sound 
source;  and  3)  the  source  term  has  been  shown  to  be  a  quadrupole  as  is  consistent  with  the  Lighthill 
formalism  that  underlies  Ribner’s  procedure. 


As  has  been  alluded  to  in  the  literature,  [12]  [8],  the  dilational  point  of  view,  despite  its  mathemat¬ 
ically  exact  derivation,  does  not  appear  to  be  practical  with  the  assumptions  that  are  required  in 
an  acoustic  analogy.  The  acoustic  analogy  formalism  requires  a  model  for  the  source  term  d.  The 
dilatation  would  be  related  to  the  solenoidal  motions;  this  is  done  by  definition  as  d  would  be  given 
in  terms  of  p^.  However,  a  convolution  must  be  taken  and  the  procedure  loses  compactness  and 
direct  relation  between  sound  source  and  a  local  (compact)  fluid  mechanical  mechanism  is  lost;  the 
integral  is  over  a  region  larger  than  that  which  contains  the  vorticity  (Howe  [6],  Ffowcs  Williams 
[12],  Crighton  [8]),  and  integrations  must  be  done  very  precisely.  This  is  consistent  with  the  scalings 
Pj  ~  t’iVj  ~  The  fact  that  this  is  not  an  easily  applicable  result  in  the  context  of  an 

acoustic  analogy  is  not  relevant;  our  purpose  is  a  more  systematic  clarification  of  the  Ribner’s  ideas 
and  a  presentation  of  consequences  of  the  dilatational  point  of  view  that  bear  relevance  to  the  DNS 
of  aeroacoustical  flows  in  fundamental  studies  such  as  [l],  [2], 


4.  Some  consequences  in  the  context  of  compressible  DNS 

The  dilatation  indicated  by  the  diagnostic  relationship  derived  above  is  the  dilatational  field  asso- 


ciated  with  the  incompressible  pressure,  p^.  The  dilatational  field  computed  in  the  DNS  of  a  sound 
generating  flow  will  contain  a  contribution  from  the  pseudo-sound  pressure,  ps,  and  a  higher  order 
contribution  from  the  acoustic  pressure,  p^.  This  can  be  shown  from  the  next  higher  order  expan¬ 
sion  of  the  continuity  equation.  In  the  source  field  the  dilatation  associated  with  the  incompressible 
pressure  (“pseudo-dilatation”)  using  the  continuity  equation,  scales  as 


1  D 


Pc* 


DV 


1  D 
PooClo  Dt 


Ps 


m 


(25) 


since  p^  ~  poo'll-  Although  is  a  measure  of  the  compressibility  of  the  flow,  it  is  not  to  be 


interpreted  as  an  acoustic  variable.  In  the  sourceless  region,  O^pa  =  0,  and  the  usual  linear  scalings 
Po.  =  PaC^  =  pc«u,’  are  relevant;  thus  w  —  ^c.  For  quadrupole  radiation  the  Lighthill  analogy,  [3], 
[8]*  [6],  gives 


, 


I^M/,  the  dilatation  associated  with  the  acoustic  field  scales  as 


w  c  Pa  «c 

Oa  ~  c-  ~  ~ - —  ' 

A  X  ^ 


(26) 


since  u^/i  ~  c/A.  Thus  d,^  ~  Mfd^.  The  acoustic  dilatation,  in  the  acoustic  regime  is,  of  course, 
related  to  the  pressure  in  the  usual  way,  da  ~  Pa  satisfies  D^da  =  0  as  does,  of  course,  d^. 
The  dimensional  set  [d^,  da]  correspond  to  the  nondimensional  set  [di,  d2].  (in  our  nomenclature  the 
subscript  1  on  the  leading  order  term  has  been  dropped). 


In  fundamental  studies  involving  the  compressible  simulation  of  aeroacoustical  flows,  [1]  or  [2], 
both  da  and  d^  are  calculated.  In  the  source  region  d^  will  dominate  and  the  dilatation  has  little 
to  do  with  the  acoustic  field  per  se;  it  is,  however,  the  source  of  the  acoustic  field.  This  is  subtly 
acknowledged  in  [l]  in  which  the  near-field  dilatation  is  not  shown  see  their  Figures  5,  6  and 
10.  The  near-field  acoustic  dilatation  could  in  principle  be  plotted  after  subtracting  dj  from  the 
computational  data.  Far  enough  away  from  the  source  region,  |a:|  >  A,  the  acoustic  portion  of  the 
dilatation  will  of  course  dominate,  da  >  d^  as  the  dilatation  associated  with  the  solenoidal  velocity 
field  through  the  pressure,  ds,  which  decays  since  d  ~  Ip^l  ~ 


5.  Conclusions 

The  intention  of  the  present  article  has  been  to  explore  and  clarify  several  issues  relating  to  Ribner’s 
[7]  dilatational  theory.  Given  the  controversy  Ribner’s  ideas  caused  some  thirty  five  years  ago,  the 
issue  appeared  to  warrant  attention  that  appears  to  be  missing  in  the  published  literature.  In 
vindication  of  Ribner’s  contention,  the  sound  source  has  been  formally  related  to  the  dilatation 
-  a  link  never  rigorously  made  in  Ribner’s  [7]  treatment.  The  analysis  has  also  indicated,  in  the 
context  of  the  flows  of  interest  in  Ribner’s  formulation,  additional  sound  source  terms. 


There  is  relevance,  in  the  present  era  of  computational  possibilities,  of  Ribner’s  ideas  as  a  diagnostic 
tool  to  be  used  to  investigate  numerical  data  from  compressible  DNS.  The  dilatational  rate  is  one 
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sound  source  in  the  near-field  and  in  the  far-field  it  is  one  of  several  adjunct  acoustic  variables 
[pa.  Pa.  da.  da]  any  one  of  which  serve  to  generate  the  linear  sound  field.  Thus  in  DNS  investiga¬ 
tions  using  dilatational  fields,  the  dilatation  rate  d.  not  the  dilatation  d.  appears  to  be  the  more 
potentially  useful  field  for  the  physical  interpretation  of  computational  results.  It  has  properties 
that  relate  to  both  the  sound  generation  and  the  sound  propagation  aspects  of  the  aeroacoustical 
problem  its  character  being  determined  by  where  it  is  evaluated. 
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Appendix:  Some  details 

The  pressure  flux  source  term,  {vkp),ku  is  now  shown  to  be  an  octupole.  The  solution  for  the 
pressure  is  ps  =  [(uiUj)  ®  Gojyij  =  The  pressure  flux  can  be  written,  using  Vk,k  =  0, 


{vkP),k  =  VkP.k  =  VkPij-Mk 


(27) 


and  the  source  term  is  seen  to  be  a  triple  divergence  and  thus  a  more  inefficient  octupole.  For 
those  not  familiar  with  this  style  of  reasoning  the  procedure  is  standard  in  aeroacoustics  a  few 
overview  comments  offered.  The  three-dimensional  Greens  function  for  the  wave  equation  is,  in  its 
far-field  approximation. 


C'oo(a;',  t';  s,  t)  =  G.x,{«  -  x  ,t  -  t')  — 


4rrc2|a;| 


S{t  —  t'  ~  \z  —  x'\c  = 


47rc^|a; 


Gs. 


(28) 


For  a  generic  source  term  S{x)  =  Siji....,ijk...  skipping  the  intermediate  steps  already  given,  the 
solution  is 


~  47rc2|®|'^*^*’-’*'^'''’''-  47rc2|a:|  ®  ^^'vk- 


(29) 


The  symmetry  in  time  and  space  of  Gg  allows  time  derivatives  to  replace  spatial  derivatives  each 
additional  spatial  derivative  contributing  a  power  of  the  fluctuating  Mach  number  as  the  symmetry 
in  {x,t)  allows  ^  Thus  Gj-,,- =  and 


Pa(®) 


_£i_  ^  ZfL  rc 


iwc^lx]  c\x\  c|a;|  c|a;| 


(30) 


Hence  the  pressure  flux  source  term,  Py.y*  is  an  order  higher  in  the  turbulent  Mach  number  than 
the  dilatation  rate. 


9 


